Energy exponents and corrections to scaling in Ising spin glasses 
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We study the probability distribution P(E) of the ground state energy E in various Ising spin 
glasses. In most models, P(E) seems to become Gaussian with a variance growing as the system's 
volume V. Exceptions include the Sherrington-Kirkpatrick model (where the variance grows more 
slowly, perhaps as the square root of the volume), and mean field diluted spin glasses having ±J 
couplings. We also find that the corrections to the extensive part of the disorder averaged energy 
grow as a power of the system size; for finite dimensional lattices, this exponent is equal, within 
numerical precision, to the domain-wall exponent 8dw- We also show how a systematic expansion 
of 9dw in powers of e~ d can be obtained for Migdal-Kadanoff lattices. Some physical arguments 
are given to rationalize our findings. 

PACS numbers: PACS Numbers : 75.50.Lk, 05.50.+q 
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I. INTRODUCTION 

Exponents and corrections to scaling play a central role 
in disordered systems. Consider for instance the depen- 
dence of thermodynamic quantities on the size L of a sys- 
tem when L — > oo. In paramagnetic systems, correlations 
are short range and thus bulk properties converge rapidly 
to their thermodynamic limits. On the contrary, in sys- 
tems such as spin glasses below T c , correlations are long 
range and so finite size corrections are "large" . In prac- 
tice, that means that even if one takes periodic boundary 
conditions, the disorder average of an intensive quantity 
converges to its thermodynamic limit slowly: finite size 
corrections go to zero as an inverse power of L. This slow 
convergence is to be contrasted with paramagnetic sys- 
tems where finite size corrections in disorder averages are 
exponentially small. One can also consider the sample-to- 
sample fluctuations of thermodynamic quantities. These 
fluctuations scale in the thermodynamic limit; an inter- 
esting open question is their limiting distribution, if any. 

In this work we consider these issues in the context of 
Ising spin glasses 0, at zero temperature and focus on 
the ground state energy. Given a distribution of disorder 
realisations (each being represented by J), the extensive 
ground state energy Ej is a random variable. Denote 
by N the number of spins in the system; N — L d for a 
d-dimensional hyper-cubic lattice. We are interested in 
the probability distribution Pn(Ej) and in how the con- 
nected moments (cumulants) of this distribution depend 
on N (or L). Denoting disorder averages by an ovcrlinc, 
we have 



Ej(L) = e L d + e x L 



(1) 



In this expression, O s is what we call the shift exponent. 
To leading order, the energy scales with the volume, while 
O s gives the leading correction to scaling. The justifica- 
tion for this notation will become clear later; in the mean 
time, we should compare to the usual notation involving 



the correction to scaling exponent ui; indeed, in the spin 
glass phase, finite size corrections of intensive quantities 
decay as a power of L: 



Ej(L)/L d = e a + ei L- 



(2) 



and so u> = d — S . 

We will also look at the width of Pm which gives the 
scaling of the instance to instance fluctuations: 



Ej(L)-Ej(L) 



1/2 



(3) 



G/ being the fluctuation exponent. Naturally, one may 
extend these types of expansions to any cumulant of P/y, 
a central question being whether P/v has a limiting shape 
as N — ► oo. One motivation for this is the fact that the 
ground state energy is an extreme statistic: there are 2 N 
configurations of spins (Si = ±1) and one is interested 
in the one of minimum energy. (Note that these 2^ ran- 
dom energies are correlated.) Little is known about such 
statistics except in a few solvable cases 0, Q . The pur- 
pose of this work is to determine Pn(Ej) numerically for 
a variety of spin glass models. From our measurements 
we extract O s , O/, etc... and compare these results to 
theory. 

The outline of this paper is as follows. We first con- 
sider expectations arising from analogies with other sys- 
tems as well as some known results. In section IIIII we 
describe the different types of spin glasses used in this 
study. Then we give in section IIVI our estimates of the 
fluctuation exponent 9 / . This is followed by a study of 
the probability distribution of the ground state energy in 
section[V] for most models we find that it becomes Gaus- 
sian in the large system size limit. Then in section IVll 
we present our results for the shift exponent S . For all 
of these studies, we compare the (i-dimensional case to 
the theoretical predictions. In section IVIII we show that 
these exponents sometimes depend on the distribution of 
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the disorder variables (the spin-spin couplings). Finally 
in section IV 1111 we discuss and conclude this work. Some 
details of the analytic computations are given in the ap- 
pendix. 



II. CLUES FROM THEORY 

What values should be expected for the two exponents 
O s and Of? Suppose we start with the random energy 
model (REM) 5| as a guide. Its ground state energy Ej 
has a Gumbel distribution in the large N limit [2j with 
a variance of 0(1). Also, the disorder averaged ground 
state energy grows as e$N + O(lniV) at large N. These 
properties show that O s = 9/ = and thus in the REM 
finite size corrections are "small" , though clearly much 
larger than in a paramagnetic system. 

Another model that can guide us is the directed poly- 
mer in a random medium (DPRM) 7]. The problem 
on a tree is very similar to the REM, although subtle 
differences appear £|. The two dimensional case (with 
one space and one "time" dimension) can also be solved 
in much detail. In these cases, perhaps surprisingly, a 
single exponent 9 describes the scaling of three a priori 
unrelated quantities: (1) the corrections to scaling of the 
disorder averaged ground state energy Ej\ (2) the size of 
sample-to-sample fluctuations of Ej; (3) the typical exci- 
tation energy of the lowest "system-size" excitation, that 
is an excitation that is macroscopically different from the 
ground state. These properties lead to the remarkable re- 
lation 9 S = 6/ = 9. The DPRM is thus described by a 
"one-parameter" scaling theory. 

It is plausible that an extension of the DPRM scal- 
ing theory may apply to Ising spin glasses. At the heart 
of such a theory, initiated by McMillan developped 
by Bray and Moore and extended by Fisher and 
Huse 10], is the exponent 9bw- This exponent is anal- 
ogous to the 9 of the DPRM, and gives the scaling of 
domain wall energies, Edw — TL 6dw ; numerical esti- 
mates give 9dw ~ 0.20 in d = 3 and 9dw ~ —0.28 in 
d = 2. (More generally, a number of other 9 exponents 
have been introduced for spin glasses; they are all asso- 
ciated with excitation energies; this is to be contrasted 
with our O exponents that are associated with ground 
state energies.) If we follow the correspondence with the 
DPRM, we expect that S = 9ow, justifying our use of 
a "theta" notation for 0^. Physically, this equality cor- 
responds to the fact that these systems are sensitive to 
boundary conditions; for some samples, these conditions 
are such that a "domain wall" must be present in the 
ground state. However the analogy with the DPRM cer- 
tainly breaks down for 0/: for any short range spin glass 
in dimension d, Wehr and Aizenman [TlT ] proved that 
0/ = d/2. This shows that the REM and the DPRM are 
not good guides for finite dimensional spin glasses. 

To have more realistic theoretical predictions, and in 
particular to preserve 0/ = d/2, it seems necessary to 
work with models having a microscopic Hamiltonian de- 



fined over configurations of N spins. One approach is to 
use hierarchical (Migdal-Kadanoff) lattices; there, ana- 
lytical computations as well as powerful numerical meth- 
ods are possible. We shall also consider mean field spin 
glasses where spins are coupled amongst one-another at 
random so there is no geometry to speak of. The P^(Ej) 
in such models can be referred to as the mean field pre- 
diction. Strangely enough, little is know about these sys- 
tems so we will have to determine their behavior numeri- 
cally. In some cases these models lead to surprizes as we 
shall soon see. 



III. MODELS AND METHODS 

We focus on three families of Ising spin glass models so 
that the effects of geometry and dimension are apparent. 
The Hamiltonian for these models is 

H = -J2 J ^iSj (4) 

<ij> 

where Si = ±1 and the Jy are quenched Gaussian ran- 
dom variables of zero mean and variance 1 (except in 
the case of the Sherrington-Kirkpatrick model which we 
discuss later). The sum X^<y> ^ s over nearest neigh- 
bor spins on a given graph having N nodes; the different 
models we consider vary simply by the nature of that 
graph. 

Our first family of models are of the Edwards- Anderson 
(EA) type 12] : the graphs are square or cubic lattices 
of linear size L (thus N = L d when the dimension is d) , 
and the edges connect nearest neighbors only; we take 
periodic boundary conditions in all directions. 

Our second family of models come from the Migdal- 
Kadanoff (MK) approach 0] where one performs a 
bond-moving real-space renormalization group. This 
procedure effectively amounts to computing quantities 
on hierarchical (MK) lattices defined by an iteration pro- 
cess (see Fig. P). The iteration takes one bond (that is 
an edge of the current graph) into b paths, each made of 
2 segments (that is edges); if r is the iteration number 
(beginning with r = 0), the "linear" lattice size L grows 
as 2 r and the volume (actually the number of edges and 
thus the number of terms contributing to the energy) 
grows as (26) r . When using I segments instead of 2 in 
each path, we have 

L = i r and N ~ {£ b) r (5) 

so that the dimension is d = In N/ In L = 1 + ln&/ln£ 
The usual choice to obtain d = 3 is I = 2 and 6 = 4, 
while £ = b = 2 corresponds to d = 2 as in Fig. ^ 

Our third family of models come from mean field, and 
here we have considered two types of graphs. First, 
we use complete graphs where all vertices are con- 
nected, corresponding to the Sherrington-Kirkpatrick 
(SK) model To have an extensive energy, one takes 
the variance of the Jy to be 1/N. Second, we also use 
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FIG. 1: Construction of a d = 2 MK hierarchical lattice. 



diluted models for which the connectivity is fixed and 
identical for all the vertices of the graph (lij. The dis- 
order ensemble then consists of the uniform distribution 
over all graphs satisfying that constraint in addition to 
the disorder ensemble in the Jy. Such an ensemble can 
be used to "model" the Euclidean case by setting the 
coordination to that of the lattice of interest. Thus to 
model the d-dimensional EA model on the hypercubic 
lattice, we set the coordination to 2d. 

For these three families of spin glass models, we will 
determine the distribution of the ground state energies. 
The hierarchical lattices allow one to write a recursion for 
Pn(Ej); because of that, there is no need to sample the 
disorder variables, we can perform the disorder average 
exactly. On the contrary, for the EA and the mean field 
cases, we must compute the actual ground state energies 
for a large number of disorder samples. For that we rely 
on a previously tested |16| algorithmic procedure where 
given enough computational ressources, the ground state 
is obtained with a very high probability for both Eu- 
clidean lattices and for random graphs as long as N is 
not too large. For our runs, we used several months of 
CPU time on Pentium III personal computers running at 
666 MHz. With this amount of CPU, we obtained high 
statistics for lattices of sizes up to L = 10 (in both d = 2 
and 3) and for mean field graphs with N up to 300. 



IV. SAMPLE TO SAMPLE FLUCTUATIONS 
AND THE EXPONENT 6/ 

A. Migdal-Kadanoff lattices 

Let us begin with the Migdal-Kadanoff approach for 
which the important prediction Qf = d/2 can be de- 
rived. To understand how this relation comes about on 
the hierarchical lattices, we construct these by "aggre- 
gation", i.e., by recursively (and hierarchically) joining 
sub-lattices together. (This procedure is to be distin- 
guished from the top-down iteration used in Fig.^) We 
work with the distribution of ground state energies E p 
(respectively E a ) subject to fixed boundary conditions, 
the spins on the ends of the lattice being forced to be par- 
allel (respectively antiparallel.) Let a 2 be the variance of 
the ground state energy at some level of the hierarchical 
construction. To go to the next level for which L will be 
£ times larger, first find the ground state energy in one 



of the b paths. Clearly, if we take the ground state con- 
figuration in each of the segments of that path, we will 
have built the ground state for unconstrained end spins. 
If the result does not give the imposed values for the end 
spins (parallel or antiparallel), one must add a "correc- 
tion" term equal to the smallest domain wall energy of 
the £ segments of that path. Thus the energy of one path 
is the sum of the ground state energies of the £ segments, 
plus one domain wall energy half of the time. Second, we 
add up the independent contributions from the b differ- 
ent paths, leading to E p and E a at the new level. The 
ground state energy is then simply min(E p , E a ). If we 
neglect the "correction" term, then a 2 at this new level is 
just £b times larger than at the previous level. At large b, 
the correction is in fact small and so it can be neglected. 
From one level to the next, the volume grows by a factor 
£b, just as <7 2 does, so the variance is linear in the lattice 
volume and thus 0/ = d/2. From our numerical study 
of these hierarchical lattices, we find that this relation 
holds also for small b and for both £ = 2 and £ = 3. (We 
did not test for larger £.) 

As mentioned in the introduction, the Wehr-Aizenman 
theorem [TT| shows that 9/ = d/2 in finite dimensional 
spin glasses. It is rather conforting that the MK approach 
also leads to this result, sustaining the belief that it is a 
useful guide for real (finite dimensional) spin glasses. 



B. Mean field models 

Next, consider the mean field prediction for Qf. Since 
there is no geometry in our mean field family of models, 
we identify L with N 1 ^ and thus 



E 2 (N)-Ej(N) 



1/2 



= a N 



Bf/d , 



(6) 



The issue here is whether mean field also predicts Qf/d — 
1/2 as it should if one believes that this approach cor- 
rectly describes the large dimensional limit of real spin 
glasses. To find out, we have performed ground state 
computations on thousands of samples of the SK model 
and of fixed connectivity spin glasses, and have extracted 
for each ensemble its associated Qf/d. 

First, consider the fixed connectivity spin glasses with 
the connectivities, z = 3,4,6 and 10. To estimate 
Qf/d, we determine when the root mean square devia- 
tion (RMS) of the ground state energy divided by N Bf l d 
becomes flat as a function of N. (Our runs were per- 
formed for 50 < N < 300.) This leads to Q f /d « 0.5 for 
z = 3, Qf/d fa 0.49 for z = 4, Qf/d « 0.48 for 2 = 6, and 
Qf/d~ 0.44 for z = 10. The drift we observe in Qf/d is 
most probably an artefact of our procedure, and simply 
corresponds to the fact that the corrections to the scaling 
in Eq. are important in our data, especially at large 
z. To get better estimates of Qf/d, we would need to 
control these corrections to scaling but our data are not 
sufficiently precise for that. Nevertheless, it seems very 
plausible that Qf/d — 0.5 in all the fixed connectivity 
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FIG. 2: Scaled root mean square deviation (RMS) of the 
ground state energy in the SK model, displayed as a func- 
tion of N~ 2/3 . Main figure: RMS divided by N 1/4 ; inset: 
RMS divided by N 1/c '. 



models. In direct analogy to what was stated previously, 
we can also say that this result sustains the belief that 
mean field models provide a useful guide to real (finite 
dimensional) spin glasses. 

Let us now continue and consider the limit of infinite 
connectivity, i.e., the SK model. Since there are 0(N 2 ) 
terms contributing to the Hamiltonian, a very simple 
minded guess would give Qf/d= 1 and thus larger fluc- 
tuations than in the other models. However the opposite 
happens, revealing that the scaling in the SK is quite 
subtle. One can get some clues from Kondor's analytic 
study |17| performed just below the critical temperature 
T c . In particular, his results have been interpreted by 
Crisanti et al. who argued that the free energy fluctu- 
ations should scale as TV 1 / 6 , which is very small compared 
to the naive estimate. A different estimate was recently 
proposed by Aspelmeier et al. who argued that free 
energy fluctuations should scale as TV 1 / 4 . Our concern 
here is the ground state energy; although there are no 
analytical calculations, it is plausible that the exponent 
for energy fluctuations at T = is the same as at that for 
free energy fluctuations at < T < T c . (Such an extrap- 
olation to T = is known to apply to the DPRM where 
the exponent of the free energy fluctuations at T > is 
equal to the one of the ground state energy fluctuations.) 
In that case, we would have Qf/d = 1/6 according to 
some authors and Qt/d= 1/4 according to others. 

What do the numerical estimates tell us about this 
question? We are aware of a study by Cabasino et al. ppl ] 
who showed beyond any reasonable doubt that Qf/d< 
0.5; in fact their best fit gives Qf/dra 0.28. When we 
consider our data, motivated by Eq. we find that the 
ratio of the left and right hand sides is compatible with 
a constant when we use 0/ ~ 0.25 and 40 < N < 150, 



but the two points at N = 200 and 300 are then below 
the others as show in Fig. [3 It is difficult to extract an 
error bar on the value of this exponent, and most likely 
the terms dropped in Eq.|S|are important just as was the 
case in the fixed connectivity models. In fact, here the 
difficulty is even more acute as the standard deviation 
increases only very slowly with N. Nevertheless, let us 
fit the standard deviation to a pure power; then using 
a x 2 analysis we obtain Qf/d = 0.24 ± 0.005 which is 
very close to the conjectured 1/4 val ue, a result further 
supported by other numerical work |2l| . Note however 
that throughout all this paper, the error bars given are 
statistical only. . . In our fit, there are 7 degrees of freedom 
and the resulting % 2 w 9.5 is not bad; nevertheless, a 
critical examination of the figure leads one to conclude 
that the actual uncertainty on Qf/d is certainly much 
greater than 0.005. Note that this estimate is a bit lower 
that that of po| and is many standard deviations away 
from the conjecture Qf/d = 1/6. However, if we impose 
Qf/d — 1/6, the data is reasonably straight as a function 
of 7V~ 2 / 3 as shown in the inset of Fig. Because of this 
"good" behavior, we cannot rule out Qf/d = 1/6. 

Given that in the dilute spin glasses we expect Qf/d — 
0.5, how does one recover the SK case as z — > oo? Begin 
by recalling that, in the SK model, the variance of the Jy 
is taken to be 1 /N to ensure that the ground state energy 
is extensive. To maintain this property in the diluted spin 
glasses, we must divide the Jy by ^fz. The fluctuations 
for large N and z in these modified spin glasses scale as 



Ej(N)-Ej(N) 



1/2 



a (z)N Q f/ d /^- 



(7) 



with Qf/d = 1/2 from what we saw previously. The 
SK limit corresponds to setting z — N, so we see that 
smooth large TV and z limits require cr (z) ~ z M where 
H is equal to the SK value of Qf/d. We saw that it was 
difficult to obtain Qf/d in the fixed connectivity models, 
but obtaining the prefactor uq(z) is even more difficult. 
Nevertheless, we have extrapolated our data for the root 
mean square deviation divided by N ^ in the different 
finite connectivity models. This leads to the estimates 
a (z = 3) w 0.67, cr (z = 4) « 0.70, a (z = 6) « 0.73, 
and Co [z — 10) ~ 0.76. This growth is very slow; it is 
compatible with the value fi — 1/6 but much less with 
the value fi — 0.25. Given the large uncertainties in our 
extrapolations however, this has to be considered as only 
very indirect evidence in favor of Qt/d = 1/6. 

We conclude that understanding the non-trivial size 
dependence of the fluctuations of the ground state energy 
in the SK model remains quite a challenge. Let us pro- 
pose here a simple argument suggesting that Qf/d = 1 /4. 
Imagine that one changes slowly all the couplings Jy = 
Jfj/yN: J°- — > J®j +5J®j, where the order of magnitude 
of SJ®j is denoted by e. If e is infinitesimal, the ground 
state remains the same. Up to what value of e will this 
be true? The change of local field induced by the change 
of couplings for a fixed configuration {S*} of the spins 
is of order V ■ SJ^SUVN — e. Using the fact that the 
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local field distribution vanishes linearly for small fields in 
the SK model, it is easy to show that the smallest local 
field is of 0(l/y/~N). Therefore, the first value of e that 
will trigger a change of the ground state is e* ~ l/\/~N. 
Flipping the spin with this low local field should lead to a 
cascade of flips that lowers the energy and thus to a new 
ground state that is "substantially" different from the 
starting one. Furthermore, the variation in the ground 
state energy when going from e = to e = e* is proba- 
bly 0(1); indeed, there are excited states of energy O(l) 
above the ground state, and one therefore expects level 
crossings to occur when e ~ 1/ \/~N. Finally, in order to 
scan the whole distribution of Jf/s, one needs e ~ 1. Us- 
ing this range, the ground state will change 1/e* ~ TV 5 
times. Since between each level crossing, the ground state 
energy randomly changes by an amount O(l), the total 
expected fluctuation of the ground state energy will be 
OiN 1 / 4 ), so that Qf/d = 1/4. More analytical work is 
obviously needed to confirm this speculative result, but 
note that at the heart of our argument lies the fact that 
the ground state of the SK model is particularly fragile: 
a relative change of order 1 / yN of the disorder is enough 
to substantially change the ground state. 



C. Edwards- Anderson models 

In the case of the finite dimensional lattices, we know 
that 8/ = d/2 holds exactly because of the Wehr- 
Aizenman theorem [llj . and recently Aspelmeier and 
Moore |22jl] found that this relation holds within replica 
theory. In spite of these theoretical results, it is instruc- 
tive to see how this equality transpires numerically. We 
thus follow the procedure used in the fixed connectivity 
models where we tested for when the rescaled RMS be- 
came size independent. In d — 3, the rescaled data show 
no obvious trend when 1.49 < 9/ < 1.60 while in d = 2 
the corresponding range is 1.00 < 9/ < 1.02. In particu- 
lar, in Fig.[3]we show these ratios when 9/ is set to d/2 
(the data displayed are from the models with Gaussian 

JijS.) 

In principle, it would be interesting to find the correc- 
tions to this scaling law. In a renormalization group pic- 
ture, corrections go as inverse powers of L. Furthermore, 
if one believes that the exponent 9 S gives the leading 
corrections to scaling for all extensive quantities, then 
those for the variance of Ej should go as L &s . To con- 
sider this possibility, we set 9/ = d/2 and then ask when 
our data for the rescaled RMS follow approximately a 
straight line when plotted as a function of with the 
expectation that lj — d — 9 S . In d = 3, we find u> > 1.5, 
to be compared with the "theoretical" value co sa 2.8 (cf. 
section lyi)! . Since 9 S is small and may in fact be zero, we 
have performed the analysis for to = 3; using the inter- 
cept of the curve with the y axis, we find that the RMS 
of Ej grows as 0.765L 3 / 2 . In d = 2, the rescaled data 
is too flat and so in practice we cannot give any sensible 
estimate of lo. But we can follow the d = 3 procedure, 
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FIG. 3: Standard deviation of the ground state energy divided 
by the square root of the volume for the EA model in d — 2 
and 3, Gaussian J^s. 



setting this time to = 2 (9 S is also relatively close to 
zero in d — 2); the corresponding fits give that that the 
RMS of Ej grows as 0.725L. The moral of this story is 
that corrections to scaling are in general very difficult to 
determine, even if the leading scaling law is known ex- 
actly. Nevertheless, it seems that the same w may very 
well describe the dominant corrections to scaling of many 
observables, as expected from the renormalization group 
picture. 

To put these last numbers in perspective, consider the 
Mattis model where the couplings are gauge transformed 
from a ferromagnet having chosen randomly on the 
positive side of a Gaussian. The total ground state energy 
of such a system is the sum of all these couplings and thus 
has the expectation value d N \J2/-k on the hypercubic 
lattice of dimension d. On the other hand its variance 
is dN[l - 2/tt]. Thus in this Mattis model, the RMS 
of the ground state energy grows as 1.04L 3 / 2 in d = 3, 
and as 0.852L in d = 2. As expected, the Mattis model 
has larger absolute fluctuations than the EA model. It is 
also appropriate to compare the relative fluctuations ay, 
that is the RMS of the ground state energy divided by 
its mean. For the Mattis model, we find a r — 0.434/L 3 / 2 
in d = 3 and a r — 0.534/i in d = 2. These should be 
compared to the values we find in the EA model: o~ r = 
0.450/L 3 / 2 in d = 3 and a r = 0.551/i in d = 2. This 
shows that the relative fluctuations are slightly smaller 
in the Mattis model than in the EA model. Although 
this is in line with what frustration should do, note that 
the size of the effect is about 3% which is very very small. 

To summarize our study of 9/ for the different spin 
glass models, we have found that all the models consid- 
ered seem to satisfy 9/ = d/2. The notable exception is 
the SK model for which we made the case that 9/ = d/A 
was likely but 9/ = d/6 was also possible; no matter 
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what, 0/ = d/2 is excluded, showing that the N — > oo 
and d — > oo limits do not commute. In all other cases, the 
variance of the ground state energy grows linearly with 
the system's volume; this is the scaling expected from a 
central limit theorem behavior when the different terms 
contributing to the ground state energy are independent. 
Our results thus tell us that these terms are only weakly 
correlated. This feature is completely missed by both the 
REM and the DPRM; although a one-parameter scaling 
picture applies to those two models, it cannot apply to 
spin glasses. Note that this central limit behavior sug- 
gests that Pn(Ej) tends to a Gaussian; we now turn to 
see whether this is the case. 



PROBABILITY DISTRIBUTION OF 
GROUND STATE ENERGIES 
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If the central limit theorem (CTL) were applicable, not 
only would we have 0/ = d/2 for the scaling of the fluc- 
tuations, but also the shape of Pn(Ej) would become 
Gaussian at large L. This behavior indeed arises for the 
MK lattices, both analytically at large b and numerically 
for all b. (In our numerical study, we find that the skew- 
ness and kurtosis Pn(Ej) decrease fast towards zero as L 
grows.) Obviously, the terms contributing to the ground 
state energy are not independent but their correlations 
are not strong enough to prevent a CLT large L scaling. 

The question we address here is whether this simple be- 
havior also holds in the other models. Let us begin with 
the mean field case. For the fixed connectivity mean field 
graphs, our data for the skewness and kurtosis of P^(Ej) 
decrease in magnitude as N increases; this decrease is 
compatible with an extrapolation to zero as N — > oo 
as illustrated in Fig. 0] (These quantities are difficult to 
measure to high precision, so this should be considered 
as only suggestive of a Gaussian limit for Pjy(Ej).) The 
SK model however is clearly in a different class: not only 
does it have 0/ ^ d/2 but also its Pm(Ej) is not Gaus- 
sian. In particular, its skewness at large N stabilizes 
around — 0.43± 0.02, while its kurtosis stabilizes around 
0.40 ±0.03. It is instructive to compare this to the values 
predicted by the REM model where P(Xj) is a Gum- 
bel distribution: there, the skewness is —1.139 while the 
kurtosis is exactly 2.4. Our estimates do not agree either 
with the values from the Fisher or Weibull universality 
classes: we thus conclude that the SK ground state en- 
ergy distribution does not belong to one of the known 
universality classes of extreme statistics [6J- 

Finally we come to the EA lattices. Given that the 
MK and diluted mean-field graphs lead to the same con- 
clusion, it will come as no surprise that our data for the 
EA lattices are also compatible with a Gaussian limiting 
shape for P^(Ej). Note that this is expected though 
not proven from the work of Wehr and Aizenman hlj, 
while it does follow from replica theory calculations |23 • 
An examination of the skewness and kurtosis of Pn(Ej) 
shows that they decay with system size, both in d = 2 



FIG. 4: Skewness of the ground state energy distribution for 
mean field and EA models with Gaussian Ji,s. 



and d — 3. Although our measurements lack precision 
when the number of spins is large, the extrapolations 
suggest zero limiting values as L — > oo as one can see 
in Fig. QJ Perhaps it is also worth noting that our data 
are quantitatively similar in the different models; for in- 
stance, for L = 8 and d = 3, we find the skewness to be 
-0.18 in EA and -0.21 in MK. 

The overall situation indicates that Pn(Ej) is Gaus- 
sian. Furthermore we checked whether the convergence 
to this Gaussian follows the central limit theorem law. 
Indeed, that law predicts for instance that the skewness 
scales as iV -1 / 2 ; thus we have plotted in Fig.0]the skew- 
ness for the different models as a function of that scale. 
The data is completely compatible with a linear conver- 
gence to zero, confirming the CLT scaling. Only in the 
case of the SK model are the variables contributing to 
the ground state energy sufficiently correlated to prevent 
a central limit theorem behavior. A physical interpreta- 
tion of this is that as soon as there is a kind of locality 
which allows one to decompose the sample into quasi- 
independent subsystems, the central limit theorem be- 
havior will appear, leading to Of = d/2 and a Gaussian 
Pn(Ej). Of course, it is not clear why this should apply 
to the diluted mean field graphs. 

VI. THE SHIFT EXPONENT 9 S 

A. Migdal-Kadanoff lattices 

We now move on and study the finite size corrections 
to the mean energy density. Following Eq. ^ the mean 
excess of the (extensive) ground state energy is expected 
to scale as L 03 . To have an idea of what this expo- 
nent should be, it is again most useful to begin with the 
hierarchical lattices. The important prediction of that 
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approach is that Q s — Odw where Odw is the domain 
wall exponent. To see why this is so, reconsider the evo- 
lution equation for the energies E p and E a as one applies 
the recursion. First, along a given path, the energy is 
the sum of the ground state energies of each of its i seg- 
ments, the sum being sometimes corrected by the energy 
of the domain wall of the weakest segment in order to 
satisfy the boundary conditions. This correction shifts 
the path's energy by 0(L 6dw ). Second, adding the en- 
ergies of the different paths does not change the scale of 
the shifts, so necessarily Q s = Odw- Naturally, we have 
confirmed this relation numerically for different values of 
I and b by performing fits to moments of P(E p , E a ) (de- 
termined with no statistical error for these hierarchical 
lattices) from which we extracted estimates for O s and 
Odw- 

When b is large, we can derive the analytical expres- 
sion for Odw on these MK lattices. Indeed, in this limit, 
P(E p ,E a ) becomes Gaussian and so one can perform a 
cumulant expansion about this Gaussian, leading to a 
1/b series. This scheme extends the work of Southern 
and Young 0| who assumed that P(E p ,E a ) was Gaus- 
sian even for finite b. For general I and b we obtain 



(8) 



When t = 2, we find a = -1.2302 and ai = -0.04573 
(see the appendix for a derivation). In Fig. [5] we plot 
the difference between the numerically obtained 0dw{^ = 
2, b) and the terms of the expansion given in Eq. [S] This 
allows us to determine numerically the next term of the 
expansion, and we find —0.045 ± 0.001/6 2 . This value 
could be obtained analytically, but we have not pushed 
the analytical calculation to that order. Note that the 
1/b expansion corresponds to an expansion in e~ d where 



d is the dimension of space. This justifies the fact that the 
1 jb expansion is quite accurate all the way down to d = 3 
(which corresponds to b = 4). Finally, when I — > oo, we 
obtain ao = 3/2. These results show that & s — Odw < 
0/, justifying the neglect of the "correction" terms in 
section HV Al from which we concluded that 0/ = d/2. 

Before going on to the mean field case, let us remark 
that the MK value for Odw is quite close the actual 
value in the EA model. If we use the standard choice 
for d — 3, I — 2 with b = 4, the MK prediction is 
Odw ~ 0.255. One can also use the choice I = 3 with 
b = 9 for which Odw ~ 0.242. These values are to be 
compared to current estimates for Odw hi the EA model, 
DW = 0.21 ± 0.02 and D w = 0.19 ± 0.02 H]. A 
similarly good comparison occurs when d = 2. 



B. Mean- field models 

We have seen that the prediction of the MK approach 
is Q s = Odw', in d — 3, this gives either 9 S = 0.25 (the 
MK value) or <d s w 0.20 (if one uses the d = 3 EA model 
values for Odw)- In the case of the mean field models, 
there is no way to introduce domain walls, and so we 
simply focus on their prediction for S . 

Consider first the finite size effects in the SK model. 
Parisi et al. [25j have computed analytically how various 
quantities converge to their large N limit. An N~ 2 ^ 3 
convergence is the general rule, though for the energy 
density they were able to compute the finite size cor- 
rection only at the critical temperature T c and on the 
de Almeida- Thouless line. Nevertheless, the natural ex- 
trapolation is that this law should apply to all T < T c , 
leading to the prediction Q s /d — 1/3. (Note that this 
prediction is very different from that of jl8j for which 
® s /d = Qf/d = 1/6.) To our knowledge, the possibil- 
ity that <d s /d = 1/3 for the SK ground state energy was 
first brought up by Palassini 26] in his numerical studies. 
We can extend his analysis with our data; identifying as 
before L with N x / d , we perform fits of the mean ground 
state energy to Eq.[T] We then find @ s /d = 0.34±0.02, in 
complete compatibility with 1/3; note that when plotted 
as a function of N~ 2 / 3 , the data is very linear start- 
ing from N = 50 (see the inset of Fig. EJ. In addi- 
tion, if we perform the fit while forcing Q s /d = 1/3, 
we find that the ground state energy density at N = oo 
is e = —0.7637 ± 0.0002, in very good agreement with 
the exact result e = -0.7633 . . . 

In contrast to what happens for &f/d, Q s /d is the 
same in the SK model and in the fixed connectivity 
models we have considered; this is illustrated in Fig. [S] 
For instance, for connectivity 6, a power law fit gives 
Q s /d = 0.35 ± 0.03 with X 2 = 9.7 for 7 degrees of free- 
dom. In fact, the value Q s /d =1/3 works very well for 
all four connectivities we studied, and we are tempted to 
consider that this value is the exact exponent. The same 
conclusion was reached by Boettcher 27] . 
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FIG. 6: Mean ground state energy density versus N~ 2/3 for 
the diluted mean field model with connectivity 6 and for the 
SK model (inset). Error bars are included. 



FIG. 7: Mean ground state energy density and the best fits 
for the d = 2 (inset) and d — 3 EA spin glass with Gaussian 
couplings. We find O s (d = 2) » -0.35 and 9 s (d = 3) as O.f. 
These value are "close" to 8dw, suggesting O s = Odw- Error 
bars are included. 



C. Edwards- Anderson models 



Now it is time to compare the MK and mean field "pre- 
dictions" to our measurements of 9 S in the EA models. 
Let us begin with the case d = 2. We fit our 4 < L < 10 
data to Eq. n where eo, e\ and <d s are ajustable param- 
eters; that gives us the estimate O s = —0.35 ±0.01. 
(The associated x 2 1S 1-9 for 5 degrees of freedom.) 
The resulting fit is displayed along with the data in 
the inset of Fig. [7| Given our statistical error and per- 
haps more importantly the systematic effects associated 
with L not being very large, this value is reasonably 
compatible with the MK prediction S = 0ow since 

e DW = -o.28 HHHi^ 

A further argument in favor of S = Odw is as follows. 
It is known that whether one modifies the boundary con- 
ditions [U or increases the size of the system [32j, the 
fractal dimension of the surface of the clusters of spins 
that change is approximately the same. It is thus likely 
that the same type of excitations are involved in deter- 
mining O s and Odw- As a consequence, we believe that 
the d = 2 EA model is described by a scaling theory 
with S = Odw, but also with 0/ ^ Odw of course. Note 
that the mean field prediction (0 S . = 2/3) is clearly off; 
however, one cannot appeal to mean field when d = 2 
because one is below the lower critical dimension. 

Let us now move on to the d = 3 EA model which 
is more challenging and has a spin glass transition at 
T c > 0. We use the same fitting function (Eq. ^) as 
before; the best fit then gives a good x 2 and a mean 
ground state energy growing as — 1.700L 3 + 1.9i 6s with 
S = 0.10 ± 0.03. This fit is displayed in the main part 
of Fig. [7] However this value of S easily changes by 0.1 
when removing some of the data points, and in fact the 
fit sometimes even leads to negative values for S . Thus 



at best we can say that S is small, somewhere between 
0.0 and 0.2. One can compare this result to "theory". 
The mean field value (0 S = 1.0) is completely ruled out, 
whereas the MK prediction S = Odw is quite acceptable 
since Odw ~ 0.2 |23l l24|. However, another possibility 
is that the discrepancy we find has a physical origin and 
that in fact S ^ Odw- Since this issue is important, 
we push the analysis a bit further as follows. Given a 
putativ e value for S , we adjust eo so that the plot of 
Ej(L) — eoL 3 versus L Qa passes through the origin. For 
S outside the range [0.05, 0.2], the data has visible cur- 
vature. In the inset of Fig. E|we show Ej(L) — eoL 3 as 
a function of i; we see that there is no clear trend, so 
even S = seems possible. Such a value could be inter- 
preted from the (mean-field-like) behavior of system-size 
excitations found in this model |33l l34| . The analogous 
analysis in the d = 2 case is shown in the main part 
of Fig. |H1 when we use the value <d s — —0.28; there the 
expected corrections to scaling work quite well. 

Of course it would be useful to have a significantly 
smaller error bar for S but we cannot go much beyond 
what we have done here: the statistical error on Ej grows 
as L d / 2 , and the amount of computation time grows still 
faster, so we cannot obtain useful information at large L. 

We have no data for d = 4, but let us briefly consider 
the published work by Hartmann 35] where he used the 
Jij = ±1 EA model. (One expects its S to be the same 
as in the Gaussian case.) Analyzing his values for the 
mean ground state energy for 2 < L < 7, we find S = 
0.2 ± 0.1; this does not compare well with his estimate 
9dw = 0.65 ± 0.04. But, if we remove the L = 2 point 
from the fit, we find S = 0.7 ± 0.2 which is in good 
agreement with Q s — Odw- We have also analyzed the 
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FIG. 9: Standard deviation of Ej divided by TV 1/4 for the z = 
3 and 6 fixed connectivity models, but with binary couplings, 
Jij = ±1. Inset: skewness of the distribution of the ground 
state energy for these same two models. 



data of Boettcher and Percus [36j , and this leads to the 
same conclusion. In summary, we cannot exclude that 
Q s < OdWi but the MK prediction 9 S = Qdw works 
surprizingly well in the finite dimensional EA models. 
On the contrary, the mean field prediction is definitely 
off, and that of the Mattis spin glass is completely wrong 
since it gives Q s = — oo. 



VII. CASE OF +/ — J COUPLINGS 

It is widely believed that exponents are universal, i.e., 
independent of the detailed microscopic nature of the dis- 
order. (Note however that there are longstanding claims 
of universality violations in spin glasses |37|.) In partic- 
ular, for the d — 3 EA, numerical computations of Odw 
confirm this to a large extent: one has Qdw = 0.21 ±0.02 
for the Gaussian [23] and Qdw = 0.19 ± 0.02 for the 
Jij = ±1 24] models. However, if T c — as arises in 
d = 2, one may expect several universality classes and 
thus some influence of the microscopic properties (i.e., 
the distribution of the Jy) on the macroscopic proper- 
ties (e.g., exponents). Note first that Qdw is known to 
be different for Gaussian and binary (Jy = ±1) cou- 
plings pE);) in the d = 2 EA model. A posteriori, that 
is not so surprising since the associated quantities (for 
instance domain wall energies) go to zero rather than to 
infinity; they can thus easily be affected by microscopic 
details. This issue can be investigated within the frame- 
work of MK lattices. We find identical values of Odw m 
the Gaussian and the J — ±1 cases whenever Qdw > 0. 
However in d — 2 and using b = £ = 2, the Gaussian 
model gives Qdw = —0.22; on the contrary, the Jij = ±1 
case leads to Odw — —00, meaning that the domain- 
wall energies decrease exponentially with the size of the 



system rather than as a power law. Thus either there 
is violation of universality (which seems unlikely to us) 
or there are several universality classes when Qdw < 0; 
Amoruso et al. [38| have given evidence in favor of the 
latter possibility since the first version of this paper was 
posted. 

Given these amendments to the scope of universality, 
the microscopic details should indeed be irrelevant for 
quantities associated with diverging energy scales. The 
surprising claim we bring forward here is that this ex- 
pectation is still too strong: a counter example is pro- 
vided by Of in the mean field fixed connectivity graphs. 
Indeed, we saw in the case of Gaussian couplings that 
Qf/d — 1/2. Now in Fig. |5|we show the rescaled stan- 
dard deviation of Ej for the case where Jy = ±1 on fixed 
connectivity graphs; the standard deviations have been 
divided by iV 1 / 4 . If Qf/d = 1/2, we should see a rapid 
divergence of the plotted values with increasing TV, but 
instead the curves are relatively flat and decreasing! In 
fact, when performing fits, we find that Qf/d is between 
1/4 and 1/5. 

To give further evidence that the Gaussian and Jij = 
±1 cases scale differently, we plot in the inset of this 
figure the values of the skewness of the distribution of Ej 
using the same models and values of N as in the main 
figure. Although our data is noisy, the skewness shows no 
sign of going to zero when N — > 00. Taken at face value, 
this means that the distribution of Ej is not Gaussian 
in the large TV limit when Jij = ±1, in sharp contrast 
to what happens in the case of Gaussian couplings. To 
drive this point home further, we find that as z increases, 
the skewness and kurtosis grow in magnitude and seem to 
approach the values we find in the SK. (But as before, our 
values become imprecise at large TV.) We thus conjecture 
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that each z is associated with a different universality class 
and that, as z — > oo, one converges to the class to which 
the SK model belongs. One way to justify this is to 
consider that in these models z is related to dimension 
rather than to lattice connectivity; the universality class 
will then change with z. To a large extent, all different 
trends give further credence to the claim that Qf/d = 1/4 
in the SK model if that is the correct value for the fixed 
connectivity models with binary couplings. 

How can one understand this "break-down" of univer- 
sality when the underlying energies diverge? Recall that 
for the Jij — ±1 spin glass model on a fixed connectivity 
graph, the local environment of a spin has no disorder 
out to finite distances: any sign of the can be gauged 
away so that all the sample-dependent fluctuations arise 
"at infinity" . On the contrary, in the Gaussian case the 
Jij fluctuations are local, leading to C^iV 1 / 2 ) fluctuations 
in the total energy. Taken at face value, this argument 
also applies to coordination fluctuations; if this is true, 
the ±1 Viana-Bray model (3^ will have 0/ = d/2. 

To summarize, the exponent O j depends on the details 
of the underlying Jij distribution even though the energy 
scale (of fluctuations) diverges. Not surprisingly, we also 
find that Pn(Ej) for that system does not become Gaus- 
sian. Finally, in spite of this major change of behavior 
when going from Gaussian Jy to binary values, we find 
that Q s /d =1/3 very precisely in both cases. 



VIII. DISCUSSION AND CONCLUSION 

Let us go over the main results of the present work 
on spin glasses. First, the sample to sample fluctuations 
of the ground state energy grow as the square-root of 
the volume in almost all models so the fluctuation ex- 
ponent satisfies 0/ = d/2. Furthermore, the distribu- 
tion over disorder of the ground state energy probably 
tends towards a Gaussian in the large volume limit as 
suggested by Wehr and Aizenman [ll|. There are two 
notable exceptions to this picture: the SK model and the 
fixed connectivity mean field models having Jij = ±1. 
In those two cases, the sample to sample fluctuations are 
much smaller, and we find 0/ ~ rf/4, even though a still 
smaller value cannot be excluded. (Note that fluctuations 
which are smaller than square root of the volume N are 
also a characteristic of the directed polymer in a ran- 
dom medium.) On the other hand, for finite dimensional 
lattices, one expects the finite density of unfrustrated re- 
gions to contribute OiN 1 / 2 ) fluctuations to the ground 
state energy. A trivial example where this is the case is 
the (unfrustrated) Mattis model; there the variance of 
the ground state energy is obviously maximal and equal 
to zN [(J 2 ) — (J) 2 ] /2 for a connectivity z. 

We also studied the exponent O s giving the corrections 
to the scaling of the average (extensive) ground state en- 
ergy. For the d = 2 and 3 EA models, we find that 
the equality O s = Odw holds within our limited preci- 
sion. This means that corrections to scaling are associ- 



ated with domain walls. However, we were not able to 
rule out S = in d = 3, leaving the door open to other 
interpretations. In the context of the MK lattices, the 
equality O s = Odw holds exactly. Physically, this equal- 
ity corresponds to the fact that in a finite fraction of the 
samples, the boundary conditions force a domain wall 
"defect" to be present in the ground state. (Note that 
the finite size corrections are indeed always positive for 
fixed and periodic boundary conditions). Furthermore, 
we were able to compute Odw analytically to order 1/6 
when the dimension becomes large. In the case of the 
mean field models (SK or diluted graphs), the situation 
is quite different: the corrections to scaling grow as N 1 / 3, . 
Since one expects Odw ~ d/2 in large dimensions, the re- 
sult Q s /d= 1/3 for these mean field models differs from 
the large dimension limit of Odw /d. 

Finally, we have exhibited examples where the expo- 
nent O / depends on the distribution of the Jij , even 
though energy fluctuations diverge with size, i.e., 0/ > 0. 

Our work suggests several paths for further studies, (a) 
Can one establish the value of Qt/d for the SK model? 

(b) Since large fluctuations of order N 1 / 2 are detrimental 
in a numerical determination of the average ground state 
energy, is there a way to substract a (computable) con- 
tribution of satisfied bonds so as to reduce the variance? 

(c) Is © s = Odw m finite dimensional spin glasses or is 
S smaller? This second possibility could follow from 
other types of excitations whose exponents are smaller 
than Odw- 
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Appendix: Exponents in the Migdal-Kadanoff 
lattices and the large dimensional limit 

In this appendix, we focus on the simplest case 1 = 2. 
To write a recursion relation, one needs to keep track of 

(r) (r) 

two energies, E p and E a that give the ground state 
energy of the MK lattice at the r-th application of the 
recursion when the two exterior spins are respectively 
parallel and antiparallel. For a 6-branch lattice with I = 
2, the ground state energies at the (r + l)-th application 
of the recursion read: 

b 

e p +1) = min 

(4 r) (!.«) + E p ] ( 2 - «). E a ] (!.«) + E { a ] (2, a)) (9) 
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E { a +1) = min 



(E^ (1, a) + £?W(2, a), fiW (1, a) + (2, a)) (10) 

where the index 1, 2 refers to the two bonds in the I di- 
rection. This equation says that along all the 6 branches, 
one has to choose the orientation of the middle spin (that 
is decimated) such that the energy contribution is mini- 
mized, given the orientations of the external spins. Note 
that all the random variables that appear in this equation 
are independent as soon as they live on different bonds. 

If one assumes that the distribution P {r \E p ,E a ) of 
energies at the r-th generation can be written for large r 
in the scaling form 

P {r) {E p , E a ) = l/a 2 f([E p - £ r ] /oy, [E a - £ r ] /a r ) 

(11) 

then it is immediate to show (using the independence of 

the branches) that £ r +i = (2b)£ r and oy+i = V%b&r, 
so that £ r ~ (2b) r and a r ~ (2b) r / 2 . Since the number 
of spins is given by (26) r , one immediately finds that the 
ground state energy is extensive and that the fluctuations 
are described by0//d=l/2. 

Define as E ( p r) - E { a r) . It is easy to show that 
obeys an autonomous recursion relation [l3[ : 

b 

A (r+i) = £ e(a) min (|A«(1, a)\, |A«(2, a)|) , (12) 



where a labels branches and 



e(a) = —sign [A(l, a)] sign [A(2, a) 



(13) 



Using the independence of the A's corresponding to dif- 
ferent branches, one finds that E r , the RMS of the dis- 
tribution of the A's, obeys the following equation: 



£ 2 +1 = 4 6£ 2 / dxg{x) / d yy 2 g(y) + x 2 / dyg(y) 
Jo Uo Jx 

(14) 

where g(.) is the distribution of A/X that is independent 
of r for large r. From this relation, one finds that S r ~ X r 



with 



A = 4b 



dxg(x) 



dyy 2 g(y) + ^ 



dyg(y) 



(15) 

where of course g(.) depends on b. The energy scale for 
flipping the relative sign of the exterior spins is £ and so 
the exponent Odw is given by 



'DW 



In A 
In 2 



(16) 



In the large dimension limit, for which b — > oo, it is clear 
that, using the central limit theorem, the distribution 
of A is Gaussian. Since <?(.) is then known, A can be 
computed from Eq. 1151 giving 



(17) 



A 2 = 0.36338 b, 
a result first obtained by Southern and Young. 

When b is large but not infinite, the first correction 
to the Gaussian is of order 1/6 (because the A have a 
symmetric distribution) and reads: 



1 



k d 4 
246 dx 1 



-*72 



(18) 



where k is the kurtosis of the initial variable, i.e., 
emin(|A(l)|, |A(2)|). To first order in 1/6, this kurto- 
sis can be computed by assuming that the A are Gaus- 
sian, and one finds k — 0.434215. Injecting the ex- 
pression of g(x) in Eq. I|15|l then gives A to order 1/6: 
A 2 = 0.36338 6-0.023035, and finally the result given in 
the main text after Eq. JSJ). 

This calculation can be extended to next order: k will 
acquire a 1/6 contribution and there will be corrections 
to g{x) of order 1/6 2 coming from the non zero sixth 
cumulant of e min (|A(1)|, |A(2)|), computed as if the A's 
were Gaussian. One could also, with more work, compute 
A for I 7^ 2. In the limit t — * oo, the problem becomes 
soluble again using the theory of extreme value statistics 
for handling the variable min(|A(l)|, |A(2)|, |A(^)|). 



[1] M. Mezard, G. Parisi, and M. A. Virasoro, Spin-Glass 
Theory and Beyond, Vol. 9 of Lecture Notes in Physics 
(World Scientific, Singapore, 1987). 

[2] Spin Glasses and Random Fields, edited by A. P. Young 
(World Scientific, Singapore, 1998). 

[3] J. Galambos, The asymptotic Theory of Extreme Or- 
der Statistics (R.E. Krieger Publishing Co., Malabar, 
Florida, 1987). 

[4] D. Carpentier and P. Le Doussal, Phys. Rev. E 63, 

026110 (2001). 
[5] B. Derrida, Phys. Rev. B 24, 2613 (1981). 
[6] J.-P. Bouchaud and M. Mezard, J. Phys. A 30, 7997 



(1997). 

[7] T. Halpin-Healy and Y. C. Zhang, Physics Reports 254, 
215 (1995). 

[8] W. L. McMillan, Phys. Rev. B 29, 4026 (1984). 

[9] A. J. Bray and M. A. Moore, in Heidelberg Collo- 
quium on Glassy Dynamics, Vol. 275 of Lecture Notes in 
Phystcs, edited by J. L. van Hemmen and I. Morgenstern 
(Springer, Berlin, 1986), pp. 121-153. 
[10] D. S. Fisher and D. A. Huse, Phys. Rev. Lett. 56, 1601 
(1986). 

[11] J. Wehr and M. Aizenman, J. Stat. Phys. 60, 287 (1990). 
[12] S. F. Edwards and P. W. Anderson, J. Phys. F 5, 965 



12 



(1975). 

B. W. Southern and A. P. Young, J. Phys. C 10, 2179 [27 
(1977). 

D. Sherrington and S. Kirkpatrick, Phys. Rev. Lett. 35, 
1792 (1975). 

C. de Dominicis and Y. Goldschmidt, J. Phys. A Lett. 
22, L775 (1989). 

J. Houdayer and O. C. Martin, Phys. Rev. E 64, 056704 
(2001), cond-mat/0105617. 
I. Kondor, J. Phys. A 16, L127 (1983). 
A. Crisanti, G. Paladin, J.-J. Sommers, and A. Vulpiani, 
J. Physique I France 2, 1325 (1992). 
T. Aspelmeier, M. A. Moore, and A. P. Young, Phys. 
Rev. Lett. 90, 127202 (2003), cond-mat/0209290. 
S. Cabasino, E. Marinari, P. Paolucci, and G. Parisi, J. 
Phys. A 21, 4201 (1988). 
M. Palassini, private communication. 
T. Aspelmeier and M. Moore, Phys. Rev. Lett. 90, 
177201 (2003), cond-mat/0211707. 

M. Palassini and A. P. Young, J. Phys. Soc. Jpn. 69, 165 
(2000), cond-mat/9910278. 
A. K. Hartmann, Phys. Rev. E 59, 84 (1999). 
G. Parisi, F. Ritort, and F. Slanina, J. Phys. A 26, 3775 
(1993). 

[26] M. Palassini, Ph.D. thesis, Scuola Normale Superiore, 



B 31, 29 (2003), cond- 



Pisa, Italy, 2000. 
S. Boettcher, Eur. Phys. J. 
mat/0208196. 
.] H. Rieger et al, J. Phys. A 29, 3939 (1996). 
<] A. K. Hartmann and A. P. Young, Phys. Rev. B 64, 

180404 (2001), cond- mat/0107308, 
i] A. C. Carter, A. J. Bray, and M. A. Moore, Phys. Rev. 

Lett 88, 077201 (2002). 
] M. Palassini and A. P. Young, Phys. Rev. Lett. 83, 5126 

(1999) , cond-mat/9906323. 
:] A. Middleton, Phys. Rev. Lett. 83, 1672 (1999). 
,] F. Krzakala and O. C. Martin, Phys. Rev. Lett. 85, 3013 

(2000) , cond- mat/0002055. 
] M. Palassini and A. P. Young, Phys. Rev. Lett. 85, 3017 

(2000) , cond-mat/0002134. 
] A. K. Hartmann, Phys. Rev. E 60, 5135 (1999), cond- 

mat/9904296. 

.] S. Boettcher and A. Percus, Phys. Rev. Lett. 86, 5211 

(2001) , cond- mat/0010337. 
] P. Mari and I. Campbell, Phys. Rev. E, 59, 2653 (1999). 
.] C. Amoruso, E. Marinari, O. C. Martin, and A. Pagnani, 

Phys. Rev. Lett. 91, 087201 (2003), cond-mat/0305042. 
[39] L. Viana and A. J. Bray, J. Phys. C 18, 3037 (1985). 



[28] 
[29] 

[30] 

[31] 

[32] 
[33] 

[34] 

[35] 

[36] 

[37] 
[38] 



